Doping quantum dimer models on the square lattice 



O 
O 

(N 

Of 

< 

OO 
(N 



CO 

c3 



I 

C 

O 

o 



> 
in 

(N 
(N 

o 

O 

c3 



c 

o 
o 



X 
S3 



Didier Poilblanc, 1, 2 Fabien Alet, 1 Federico Becca, 1,3 Arnaud Ralko, 2 Fabien Trousselet, 1 and Frederic Mila 2 

1 Laboratoire de Physique Theorique, CNRS & Universite de Toulouse, F-31062 Toulouse, France 
2 Institute of Theoretical Physics, Ecole Polytechmque Federate de Lausanne, BSP 720, CH-1015 Lausanne, Switzerland 
CNR-INFM Democritos National Simulation Centre and International 
School for Advanced Studies (SISSA), 1-34014 Trieste, Italy 
(Dated: February 3, 2008) 

A family of models is proposed to describe the motion of holes in a fluctuating quantum dimer 
background on the square lattice. Following Castelnovo et al. [Ann. Phys. (NY) 318, 316 (2005)], 
a generalized Rokhsar-Kivelson Hamiltonian at finite doping which can be mapped on a doped 
interacting classical dimer model is constructed. A simple physical extension of this model is also 
considered. Using numerical computations and simple considerations based on the above exact 
mapping, we determine the phase diagram of the model showing a number of quantum phases 
typical of a doped Mott insulator. The two-hole correlation function generically exhibits short- 
range or long-range algebraic correlations in the solid (columnar) and liquid (critical) phases of 
the model, respectively. Evidence for an extended region of a doped VBS phase exhibiting holon 
pairing but no phase separation is given. In contrast, we show that hole deconfmement occurs in 
the staggered dimer phase. 

PACS numbers: 75.10.-b, 75.10.Jm, 75.40.Mg 



Soon after the discovery of cuprate superconductors 
with high critical temperatures, Anderson suggested that 
the Resonating Valence Bond (RVB) state is the rele- 
vant insulating parent state that becomes superconduct- 
ing under (arbitrary small) hole doping yj. Such a state 
can alternatively be viewed as a spin liquid (SL), as it 
has no magnetic order and it does not break any lattice 
symmetry. Since then, the search for exotic SL in micro- 
scopic or effective models has been very active. 

In quantum spin models, where magnetic frustration 
suppresses long-range magnetic order, spin liquids of- 
ten compete with quantum disordered states named "Va- 
lence bond solids" (VBS) which break translation sym- 
metry 0- This is e.g. the case in the frustrated Heisen- 
berg model with extended-range antiferromagnetic (AF) 
interactions In a VBS, nearest-neighbor spins pair 
up in bond singlets which order e.g. along columns or 
in a staggered arrangement. Hole doping has also been 
extensively studied in Mott insulators Q and AF fluctua- 
tions have been identified as the glue for pairing. Uncon- 
ventional pairing upon doping models exhibiting a VBS 
ground state has also been found [5j. 

In a pioneering work Rokhsar and Kivelson introduced 
a quantum dimer model (QDM) , a Hamiltonian acting in 
the space of two-dimensional fully packed dimer configu- 
rations Q (generically called |c)). The dimer interaction 
V and the dimer-flip process J are schematically depicted 
in Fig.^a). As discussed in Ref.|3, the QDM can be con- 
sidered as the simplest effective model to describe quan- 
tum disordered phases similar to the pseudo-gap phase of 
the cuprate superconductors. In that respect, one might 
think of the dimer flip term as originated directly from 
super-exchange between copper spins. At the special 
point V/J = l, named Rokhsar-Kivelson (RK) point, the 
ground state (GS) is exactly known and can be mapped 
onto the partition function of a classical dimer model . 
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FIG. 1: (Color online) Pictures of the two quantum processes 
considered in this paper, (a) Dimer flip J within a plaquette. 
A dimer repulsion V is defined for all flippable plaquettes. 
(b) Hole hopping along a plaquette diagonal. In this process 
a dimer "rotates" from a vertical (horizontal) to a horizontal 
(vertical) bond. The hole-dimer repulsion Vhd is defined on 
all the dashed lines. 



On the square lattice, the dimer-dimer correlations are 
algebraic, decaying as 1/r 2 . This "algebraic SL" at the 
RK point is believed to be rather singular on the VJJ 
axis since, as shown e.g. by numerical calculations Q, 
the GS is a VBS on both sides of it, a staggered phase for 
V/ J > 1, a columnar phase at attractive V (i.e. V < 0), 
separated from the RK point by a small region of plaque- 
tte phase. The case of non-bipartite lattices, where the 
RK point has a gapped GS that shows fractional excita- 
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FIG. 2: (Color online) Mean-squared separation Rhh between 
two holes in the t-J-V model as a function of the ratio V/J. 
Calculations are done on a periodic 6x6 cluster for two values 
of the hopping t/J. Inset: Size dependence of Rhh obtained 
from GFMC onixl clusters for t/J = 0.5 and two typical 
values of V/J showing hole-hole deconfinement and confine- 
ment, respectively. 



tions 0, is also of great interest. Doping was introduced 
in Ref. |6| and studied further later by Syljuasen [9J who 
computed dimer correlations and the energy of two static 
monomers in a background of dimers and dynamic holes. 
However, hole correlations of the dynamic holes them- 
selves, have not been investigated so far. 

In this article we construct various models of doped 
quantum dimers with the aim to study these correlations 
as accurately as possible. Like the more "microscopic" 
t-J model, the lightly doped QDM also provides a real- 
istic description of relevant quantum disordered phases 
(which, in fact, would be stabilized only at finite dop- 
ing in the t-J model), while being much easier to handle 
numerically jlfjj . A generalized RK Hamiltonian is in- 
troduced at finite doping. This model can be mapped 
onto a doped interacting classical dimer model enabling 
an efficient use of classical Monte Carlo (MC) ^l). In 
addition, it offers a controllable parameter (the effective 
temperature) to smoothly tune the system from a VBS 
to a liquid phase (even at zero doping). We also extend 
this model to an enlarged physical space away from this 
so-called RK axis, where such a mapping is no longer 
valid 12] and where full quantum computations such as 
Lanczos Exact Diagonalisation (ED) and Green's func- 
tion Monte Carlo (GFMC) are required. We provide ev- 
idence for hole deconfinement in the algebraic dimer 
phase. However, we also argue that phase separation 
occurs for low hole kinetic energy and, lastly, provide a 
complete phase diagram of the model. 

Let us first assume that holes are introduced by pairs 
on some of the dimer bonds. Next, the simplest way 
to account for their motion 0, is to consider processes 
like the one depicted in Fig.^b) where a hole hops along 
a plaquette diagonal with some amplitude t. Note that 
here holes are really thought of as new charge degrees of 



freedom originating e.g. from doping a Mott insulator. A 
general form of doped QDM which operates in a Hilbcrt 
space with a fixed hole number can then be written as 

ff = $>|c>(c|-j£ \c')(c\-t ]T \c")(c\, (1) 



(c,c') 



(c,c") 



where the sum over |c) includes all configurations with ar- 
bitrary hole positions. The sum over |c) and |c') extends 
on all pairs of doped dimer coverings (only) differing by a 
single plaquette-flip [shown on Fig. ^a)]. Similarly, the 
sum over |c) and |c") extends to all pairs of doped dimer 
coverings differing (only) by a single hole hopping along 
a diagonal of a plaquette and a single dimer "hop" from a 
vertical (horizontal) bond to a horizontal (vertical) bond. 
Such a process between two configurations |c) and |c") 
is depicted in Fig. ^b). As shown below, the role of the 
diagonal energies e c is crucial and various choices will be 
discussed. Whether mobile holes remain confined [l4| or 
not (or whether they form bound states) is the central 
issue of this study. 

Following Ref. [y, we first start with the simple dimer 
interaction introduced above, namely e c — e° c — VN C 
where N c corresponds to the number of flippable plaque- 
ttes in configuration | c) . The properties of two holes [l5j 
in such a t-J-V model are studied here by ED and 
GFMC [lg on finite clusters and the mean-squared hole- 
hole distance i?hh = \l (r 2 ) is shown in Fig.|21as a func- 
tion of the dimer repulsion V/J for fixed ratios t/J. A 
very abrupt variation is observed at V/ J = 1 (especially 
at small t) showing a clear confinement [l4| at V/J < 1, 
as expected in a VBS, and a large value of the hole- hole 
separation for V/J > 1, which scales linearly with sys- 
tem size (see inset of Fig|2J). In this case, the undoped 
system is also a VBS, but a pure staggered dimer state 
with no quantum fluctuations. It is therefore easy to see 
that the two holes can freely move away from each other 
in opposite directions along the same diagonal creating 
a string of dimers at 90° from the background at no en- 
ergy cost 0] . On the contrary, in a plaquette or colum- 
nar VBS the energy cost grows linearly with the string 
length leading to confinement. Note that it was argued 
that confinement is lost for a sufficiently high fraction of 
holes (monomers) [jj. 

In order to construct more general doped QDM, we 
briefly re-examine the undoped case (i.e. zero doping for 
which the t term of Eq. Q is irrelevant) and introduce a 
simple construction that extends the undoped RK point 
to an infinite axis. Following Castelnovo et al. we 
define the J-/3 Hamiltonian by introducing new diagonal 
energies as, 

e c = ef? = V Y, exp {~(3V ci (N c , ~ N c )} (2) 

c'(c) 

where the sum extends on the dimer covering |c') dif- 
fering from \c) by a single plaquette-flip [shown on 
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FIG. 3: Schematic picture of the parameter space of the 
models considered here; (a) J-0 model (undoped) and (b) t- 
J-/3 model (doped). The RK axes (thick lines parametrized 
by B) are defined by J = 1 (undoped case) and by J = t = 
Hid = 1 (doped case) . The points P and R correspond to the 
origins at /3 = on these axes. 



Fig. ma)]. I Q tne ~* limit, the expression of e^ lp 
reduces to and the significance of V becomes clear. 
Hereafter, V = 1 sets the energy scale. V c i corresponds 
to a classical dimer-dimer interaction. We restrict our- 
selves to the attractive case and use units for 8 such 
that V c \ = — 1. A cartoon of the Hamiltonian mani- 
fold parametrized by (J, 8) is shown in Fig. Eta). For 
,/ = 1 it is easy to check that the GS is simply given by 
^Ec ex P(-|/5K;l-/Vc)|c) with energy E = 0, where 
the normalization factor Z = ^2 C exp (— 8V C \N C ) can 
be considered as a partition function of a classical in- 
teracting dimer model 0. We have checked by ED 
data of a 8 x 8 cluster within its fully symmetric space- 
group irreducible representation, that the specific heat 
[defined as 8 2 {{N 2 ) - (iV c ) 2 )] is very close to the MC 
results obtained for a very large cluster (not shown). 
The model displays a Kosterlitz-Thouless (KT) transi- 
tion [2fJ at p = Pkt — 1.536 between a critical phase 
at 8 < (3kt (with 8 varying exponents) and a colum- 
nar dimer phase |l9j . Independently from our investiga- 
tion, a similar mapping was derived and the properties of 
this critical phase were investigated with transfer-matrix 
techniques |21j . 

To investigate the expected confincmcnt- 
deconfinement transition 0] f° r holes at the KT 
transition, let us now generalize the construction by 
Castelnovo et al. to finite doping. We define the 
t- J-/3 Hamiltonian by adding a second diagonal term to 



the one J3J of the 3-8 Hamiltonian, e c = e 



- .flip 



hop 



£c° P = Vm ]T exp{-~3V cl (N c „ - N c )} 

c"(c) 



(3) 



where the sum now extends on the doped dimer cov- 
erings \c ) connected to c) by a i-process [see again 
Fig. Gib)]. Note that this new term, as the i-term of 
JTJ, scales like the hole concentration. Vhd is a new en- 
ergy scale which naturally makes sense in the 8 — > limit 
discussed below. The term is of central importance 
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FIG. 4: (Color online) Mean-squared separation between two 
holes in the t-J-/3 model as a function of 0. ED data obtained 
on a 6 x 6 periodic cluster are shown for several values of the 
hopping t and J = 1. The dashed line denotes the location 
of the KT transition at (5 KT ~ 1.536. For t = 0.5 and t = 0.6 
the "kinks" signal the appearance of a two-hole bound state. 



as for t = Vhd and J = 1 the GS can again be written as 
-jm J2c ex P ( — \@Vc\Nc)\c) ■ In analogy with the undoped 
case, we can again define a RK axis which runs along an 
orthogonal direction to the three-dimensional parameter 
space spanned by J, t, and Vhd- This RK axis is given 
by J = 1, t = Vhd and parametrized by 8. A cartoon 
picture of this set of Hamiltonians is shown in Fig. E^b) . 
Since the hopping term (for t ^ 0) couples all topological 
symmetry sectors the GS with energy E = becomes 
unique. Interestingly, the procedure followed here can 
be generalized to more complicated dimer or hole kinetic 
off-diagonal processes. 

The 8 — ► limit is of special interest. As seen above, 
the diagonal term (J2J reduces to a dimer-dimer repulsion 
of magnitude V (set to 1) and the undoped RK point is 
recovered for J = 1 (when holes are not present). Sim- 
ilarly, the second diagonal contribution J2J reduces to 
a dimer-hole interaction on a plaquette [with definition 
given pictorially in Fig. db)] of magnitude Vhd- The 
8 = limit therefore gives rise to a large class of physical 
Hamiltonians parametrized by arbitrary magnitudes of 
J, t and Vhd (measured in units of V = 1). A complete 
investigation of this model is left for a future study 
and, in the following, we restrict ourselves to J = V (= 1 
for convenience) so that the "distance" from the RK-axis 
will be controlled by the deviation of t from Vhd- 

We first start with the case of two holes in the t-3-8 
model. The hole-hole correlations have been computed 
by ED of a 6 x 6 cluster for arbitrary t and, for conve- 
nience, for Vhd — 1 ■ Results are shown in Fig. 0] The 
ED results show a rather smooth variation of i?hh across 
the transition at Pkt- However, the finite size scaling 
at t = 1 obtained by classical MC in Fig. [5{a) shows a 
clear qualitative change of behavior at the KT transi- 
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FIG. 5: (Color online) Results for 2 holes in the t-J-/3 model 
(J = V(— 1)) for Vhd = V(= 1). Size dependence of Rhh 
by (a) classical MC at t = 1 for various values of (3 (log-log 
scale), (b) GFMC at j3 = 0.8 and t = 0.5 (linear scale), (c) 
Conjectured phase diagram for two holes in the thermody- 
namic limit for J = V = Vhd(~ 1) vs (3 and t. The dashed 
region corresponds to the confined phase for f3 > (3kt- 



critical 



| V hd /t doped VBS ^ V hd /t 



doped SL 

(deconfined) 




(a) 




FIG. 6: Phase diagrams of the t-J-/3 for V = J vs hole 
density and Vhd/t ratio . "PS" stands for "phase separation", 
(a) Phase diagram for (3 < /3kt; (b) Phase diagram for f3 > 
/3kt- The boundary between the confined (doped VBS) and 
deconfined regions is schematic. 



tion: while Rhh remains finite in the confined phase, it 
diverges as a power law in the critical phase in agree- 
ment with Ref. [p| For t = 0.5 and t = 0.6 we observe 
a kink in the ED data of Fig. HI Moreover, for t < 0.25 
i?hh remains always very close to 1, even when (3 — > 
(not shown). This signals the appearance of a two-hole 
bound state within the critical SL phase for (3 below /3kt- 
This scenario is supported by the GFMC data on larger 
systems shown in Fig. IJb) . We have checked by GFMC 
that this bound state, in fact, persists up to t = 1 (=Vhd)- 
These results are summarized in a phase diagram for two 
holes in Fig.|!Jc). 

Lastly, using simple arguments, we construst the phase 



diagram of the doped QDM for V — J (and bosonic 
holes). The previous finding that two holes pair up in 
some region of parameter space implies that either (i) 
a Q=2e superconducting state or (ii) a phase separated 
state occurs at finite doping. For Vht — t, the inverse 
compressibility k _1 (proportional to the second deriva- 
tive of the energy per site w.r.t the hole density) iden- 
tically vanishes since the GS energy vanishes on the RK 
axis for all number of holes. In addition, since Vhd acts 
effectively as an attraction between holes (or dimers), 
k" 1 should be a monotonous function of the ratio Vhd/t 
as checked numerically ■ These simple considerations 
imply that KT 1 < (> 0) for t < Vhd {t > V hd ) and, 
then, the phase separation boundary is exactly given by 
t = Vhd for all hole density and (3. Fig.[fJ]shows the phase 
diagram of the model for J = V as a function of t, Vhd 
and hole density. For [3 < Pkt, holons are deconfined in 
the SL phase. We also expect a similar phase diagram 
for (3 > /3kt (Fig. EJb)) although the critical line (at 
iih = 0) should be replaced by an extended (hole pair) 
region with VBS columnar order which survives up to 
a critical doping. However, the boundary of the phase 
separated region remains unchanged. The exact curve 
for the VBS/SF boundary remains to be investigated in 
more details. 

In conclusion, we have introduced a class of simple 
doped QDM on the square lattice which, we believe, pro- 
vide insights on the physical quantum disordered phases 
of real materials. An exact mapping onto a doped classi- 
cal dimer model (characterized by an inverse temperature 
(3) can be realized along a one-dimensional manifold of 
this multi-dimensional space. Using numerical computa- 
tions and simple considerations based on the above exact 
mapping, we determine the complete phase diagram of 
the model (at V = J) showing a number of interesting 
physical phenomena that, we believe, could be generic in 
the vicinity of a Mott insulator, beyond the framework 
of QDM. It is found that two doped holons are confined 
in the columnar VBS phase. At finite doping, evidence 
for an extended region of a metallic VBS phase exhibit- 
ing holon pairing without phase separation is given. In 
contrast, we find that the algebraic dimer-dimer correla- 
tions of the (undoped) critical (quasi-ordered) phase do 
not provide the confinement of two injected holons. Fur- 
thermore, as shown in the case of the staggered dimer 
phase, we point out that holon deconfincmcnt can even 
occur in a VBS. 
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